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Abstract 


Two  statistical  techniques  arc  presented  for  comparing  temporal  pulse  trains  received  from 
an  Electronic  Support  receiver  to  a  priori  templates  representing  the  expected  behaviour 
of  radar  signals.  The  techniques  implement  a  specialized  cross-correlation  procedure  and 
incorporate  the  effects  of  non-cumulative  jitter.  A  theoretical  basis  is  presented  in  which 
the  signal  of  interest  is  modeled  to  be  present  in  a  noisy  channel.  The  channel  noise  may 
include  the  failure  to  detect  some  expected  pulses  along  with  the  detection  of  unexpected 
pulses.  The  methods  arc  compared  to  the  more  straightforward  approaches  based  on  cross¬ 
correlation  histograms  and  a  simulated  example  scenario  is  presented. 


Resume 


Le  document  presente  deux  techniques  statistiques  permettant  de  comparer  des  trains  d'im- 
pulsions  temporelles  obtenus  a  parti r  d’un  recepteur  de  soutien  electronique  (SE),  avec  des 
modeles  a  priori  representant  le  comportement  attendu  de  signaux  radar.  Ces  techniques 
comprennent  une  procedure  de  correlation  croisee  specialisee  et  integrent  les  effets  de  gigue 
non  cumulatifs.  Le  document  presente  une  base  theorique  pour  la  modelisation  du  signal 
d’interet  dans  un  canal  bruyant.  Le  bruit  du  canal  peut  etre  du  entre  autres  a  la  non-detection 
d'impulsions  attendues  aussi  bien  qu’a  la  detection  d'impulsions  inattendues.  Les  methodes 
sont  comparees  aux  approches  plus  simples  fondees  sur  des  histogrammes  de  correlations 
croisees,  et  un  scenario  simule  est  presente  a  titre  d’exemple. 
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Executive  summary 


The  field  of  “Electronic  Support”  (ES)  depends  heavily  on  the  availability  of  a  robust  mech¬ 
anism  with  which  to  analyze  sequences  of  pulses.  In  general,  the  complete  stream  of  pulses 
received  by  an  ES  system  is  composed  of  multiple  signals  exhibiting  diverse  features,  both 
periodic  and  aperiodic,  that  arc  observed  through  a  common,  imperfect  channel.  The  com¬ 
plexity  of  the  complete  ES  problem  has  lead  to  the  development  of  numerous  techniques 
designed  to  separate  and  identify  the  constituent  signals. 

Pulses  trains  are  typically  analyzed  and  matched  against  some  form  of  a  priori  informa¬ 
tion  in  order  to  determine  the  identity  and/or  behaviour  of  one  or  more  uncooperative  radar 
signals.  The  techniques  for  temporal  pulse  train  analysis  fall  primarily  into  two  concep¬ 
tually  distinct  categories.  The  first  and  perhaps  the  more  traditional  technique  is  to  base 
the  analysis  on  generic  periodicities  that  occur  in  the  pulse  train.  Autocorrelation  methods 
arc  commonly  used  to  identify  periodic  signals  and  arc  often  implemented  using  so-called 
“dclta-t  histograms”.  Such  techniques  lead  ultimately  to  an  estimation  of  the  parameters 
associated  with  identified  pulse  patterns  and  a  comparison  of  those  estimates  with  some 
a  priori  look-up  table.  It  has  been  observed  that  these  techniques  exhibit  numerous  short¬ 
comings;  amongst  them  is  the  inability  to  detect  complex  pulse  sequences  and  patterns  with 
very  long  periods  (such  as  //-level  staggers,  where  N  1). 

A  second  approach  to  pulse-train  analysis  is  to  search  directly  for  specific  sequences  of 
pulses  using  a  priori  patterns  and  parameters.  This  approach  has  the  advantage  of  being 
somewhat  more  flexible  than  the  autocorrelation  approach  since  the  patterns  of  interest  can 
be  quite  complicated,  and  need  not  be  periodic.  This  approach  is  typically  implemented  us¬ 
ing  some  type  of  cross-correlation  between  the  observed  pulse  sequence  and  the  templates 
in  order  to  detect  a  good  match. 

In  this  memo,  two  alternative  cross-correlation  techniques  arc  presented  based  on  proba¬ 
bilistic  signal  models.  It  is  argued  that  these  approaches  arc  theoretically  well-suited  to 
account  for  such  effects  as  the  uncertainties  in  pulse  arrival  times  (whether  induced  by 
deliberate  jitter  or  by  receiver  and  channel  imperfections)  and  other  channel  impairments 
such  as  the  effects  of  missing  and  spurious  pulses.  The  methods  arc  compared  to  a  cross¬ 
correlation  histogram  using  a  simulated  pulse  scenario  and  the  results  are  reported. 
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Sommaire 


En  matiere  de  soutien  electronique  (SE),  il  est  tres  important  de  disposer  d’un  mecanisme 
robuste  d’analyse  de  sequences  d’impulsions.  En  general,  le  train  d'impulsions  complet 
rccu  par  un  systeme  SE  se  compose  de  multiples  signaux  presentant  diverses  caracteristiques, 
tant  periodiques  qu’aperiodiques,  qui  sont  observees  dans  un  canal  comrnun  imparfait.  La 
complexite  de  E  ensemble  du  probleme  SE  a  conduit  au  developpement  de  nombreuses 
techniques  en  vue  d’isoler  et  d’identilier  les  signaux  constitutifs. 

Les  trains  d'impulsions  sont  typiquement  analyses  et  mis  en  correspondance  avec  une  cer- 
taine  forme  d’information  a  priori  afin  de  determiner  l'identite  et/ou  le  comportement  d’un 
ou  plusieurs  signaux  radar  non  cooperatifs.  Les  techniques  d’analyse  de  trains  d’impul¬ 
sions  temporelles  se  divisent  en  deux  categories  distinctes  du  point  de  vue  conceptuel. 
Les  techniques  de  la  premiere  categorie,  qui  sont  peut-etre  les  plus  traditionnelles,  font 
reposer  l’analyse  sur  des  periodicites  generiques  qui  se  produisent  dans  le  train  d’impul¬ 
sions.  Des  methodes  d’ autocorrelation  sont  couramment  utilisees  pour  identifier  les  signaux 
periodiques  et  sont  souvent  mises  en  oeuvre  a  l'aide  d'histogrammes  dits  <C  delta-tau  3>. 
Ces  techniques  permettent  finalement  d’estimer  les  parametres  associes  aux  configurations 
d'impulsions  et  de  comparer  ces  estimations  avec  les  donnees  d’une  table  de  recherche  a 
priori.  On  a  observe  qu’elles  presentent  de  nombreuses  lacunes  ;  entre  autres,  elles  ne  per¬ 
mettent  pas  de  detecter  des  sequences  d'impulsions  complexes  et  des  configurations  a  tres 
longues  periodes  (p.  ex.  decalages  a  N  niveaux,  ouN  1). 

Une  deuxieme  facon  d’analyser  les  trains  d'impulsions  consiste  a  rechercher  directement 
des  sequences  d’impulsions  specifiques  en  utilisant  des  configurations  et  des  parametres 
a  priori.  Cette  approche  a  l’avantage  d’etre  un  peu  plus  souple  que  les  techniques  d’auto- 
correlation,  car  les  configurations  d’interet  peuvent  etre  tres  compliquees,  et  elles  ne  sont 
pas  necessairement  periodiques.  Elle  fait  en  general  appel  a  un  certain  type  de  correlation 
croisee  entre  la  sequence  d’impulsions  observee  et  les  modeles  pour  la  detection  de  corres- 
pondances  pertinentes. 

Le  document  presente  deux  techniques  de  correlation  croisee  fondees  sur  un  modele  de 
signal  probabiliste.  On  soutient  que  ces  approches  offrent  des  bases  theoriques  satisfaisantes 
pour  rendre  compte  d’effets  tels  que  les  incertitudes  dans  les  temps  d’arrivee  des  impulsions 
(dues  soit  a  une  gigue  deliberee  ou  a  des  imperfections  du  recepteur  et  du  canal)  et  d’ autres 
degradations  du  canal,  comme  les  effets  d'impulsions  manquantes  ou  non  essentielles.  Les 
methodes  sont  comparees  a  un  histogramme  de  correlations  croisees,  a  l’aide  d’un  scenario 
d'impulsions  simule,  et  les  resultats  obtenus  sont  presentes. 
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1  Introduction 


The  field  of  “Electronic  Support”  (ES)  depends  heavily  on  the  availability  of  a  robust  mech¬ 
anism  with  which  to  analyze  sequences  of  pulses.  In  general,  the  complete  stream  of  pulses 
received  by  an  ES  system  is  composed  of  multiple  signals  exhibiting  diverse  features,  both 
periodic  and  aperiodic,  that  arc  observed  through  a  common,  imperfect  channel.  The  com¬ 
plexity  of  the  complete  ES  problem  has  lead  to  the  development  of  numerous  techniques 
designed  to  separate  and  identify  the  constituent  signals. 

Pulse  trains  arc  typically  analyzed  and  matched  against  some  form  of  a  priori  information 
in  order  to  determine  the  identity  and/or  behaviour  of  one  or  more  uncooperative  radar  sig¬ 
nals.  The  techniques  for  temporal  pulse  train  analysis  fall  primarily  into  two  conceptually 
distinct  categories.  The  first  and  perhaps  the  more  traditional  technique  is  a  “data-driven” 
technique  whereby  the  analysis  is  based  on  generic  periodicities  that  occur  in  the  pulse 
train.  Autocorrelation  methods  arc  commonly  used  to  identify  periodic  signals  and  arc 
often  implemented  using  so-called  “dclta-x  histograms”  [1,  2,  3],  Such  techniques  lead 
ultimately  to  an  estimation  of  the  parameters  associated  with  the  identified  pulse  patterns 
and  a  comparison  of  those  estimates  with  some  a  priori  look-up  table.  As  discussed  in  [4], 
these  techniques  exhibit  numerous  shortcomings,  amongst  them  is  the  inability  to  detect 
complex  pulse  sequences  and  patterns  with  very  long  periods  (such  as  /V-lcvcl  staggers, 
where  /V>1). 

A  second  approach  to  pulse-train  analysis  is  known  as  the  “model-driven”  approach  whereby 
one  searches  directly  for  specific  sequences  of  pulses  that  arc  known  to  be  associated  with 
particular  systems,  using  a  priori  patterns  and  parameters.  This  approach  has  the  advan¬ 
tage  of  being  somewhat  more  flexible  than  the  autocorrelation  approach  since  the  patterns 
of  interest  can  be  quite  complicated,  and  need  not  be  periodic.  The  most  straightforward 
implementation  of  such  an  approach  is  to  perform  some  form  of  cross-correlation  between 
the  observed  pulse  sequence  and  the  templates  in  order  to  detect  a  good  match.  An  ex¬ 
ample  of  this  approach  is  that  of  Elton  [5]  who  advocates  a  “cross-correlation  histogram”. 
Other  related  techniques  have  been  explored  as  paid  of  a  collaboration  under  The  Technical 
Cooperation  Panel  [6], 

In  this  memo,  two  alternative  cross-correlation  techniques  are  presented  based  on  proba¬ 
bilistic  signal  models.  It  is  argued  that  these  approaches  are  theoretically  well-suited  to 
account  for  such  effects  as  the  uncertainties  in  pulse  arrival  times  (whether  induced  by 
deliberate  jitter  or  by  receiver  and  channel  imperfections)  and  other  channel  impairments 
such  as  the  effects  of  missing  and  spurious  pulses.  The  methods  are  compared  to  the  cross¬ 
correlation  technique  of  [5]  using  a  simulated  pulse  scenario  and  the  results  arc  reported. 

2  Gaussian  scoring  profiles 


In  this  section,  a  technique  is  presented  to  compute  a  probabilistic  measure  of  consistency 
between  an  observed  pulse  sequence  and  an  a  priori  signal  template.  The  signal  of  interest 


DRDC  Ottawa  TM  2004-220 


1 


is  represented  using  a  Gaussian  statistical  model  with  stationary  channel  impairments. 


2.1  Stationary  signal  model 

A  Gaussian  transmitted  pulse  sequence  model  is  proposed.  The  model  consists  of  two 
components:  a  “pulse  template”  T  and  a  time-of-arrival  variance  a2.  The  pulse  template  is 
a  set  of  nominal  times,  denoted  by  T  =  {t\ ,  to,  ■  ■  •  ,  hv}  at  which  a  pulse  may  be  transmitted, 
modulo  some  absolute  time  offset.  The  time-of-arrival  variance  a2  indicates  the  amount  of 
non-cumulative  jitter  that  may  be  associated  with  each  pulse  in  T .  Each  pulse  time  in  T  is 
known  as  a  “template  element”  and  the  number  of  such  elements  is  denoted  by  |  T  |  =  N. 

A  simple  statistical  model  of  the  transmitted  pulse  train  is  given  as  follows.  If  a  pulse 
sequence  “T  starts  at  time  offset  t”  then  the  transmitted  pulses  form  a  set  of  normally 
distributed  random  variables  denoted  by  X  =  {ij  ,*2,  •  •  •  , Tv } ■  A  model  of  “non-cumulative 
jitter”  [7]  is  adopted  so  that  the  transmitted  pulse  times  arc  independent  random  variables 
with  Gaussian  distributions 

xj~9t{tj  +  x,a2).  (1) 

The  mean  values  of  these  distributions  are  modeled  to  be  £ (xj)  =tj  +  T  and  the  variances 
arc  “E(xj  —  (tj  +  x) ) 2  =  a2.  For  simplicity,  we  suppose  that  all  of  the  transmitted  pulses  arc 
contained  within  some  observation  time  interval  whose  total  time  duration  is  denoted  by  D. 

The  transmitted  pulse  set  X  is  randomly  mapped  into  a  set  of  received  pulse  times  A  by  a 
stationary  channel  impairment  process.  Each  pulse  in  X,  independently  of  the  other  pulses, 
has  a  probability  P(\  of  being  detected  and  appearing  in  the  set  X ;  conversely,  the  probability 
of  the  pulse  being  missed,  and  not  appearing  in  X,  is  1  —  f\\ .  In  addition,  the  detected  pulses 
arc  interleaved  with  spurious  pulses.  For  simplicity,  suppose  that  the  spurious  pulses  arc 
distributed  according  to  a  Poisson  process  with  pulse  density  ps  over  a  the  observation 
interval.  The  resulting  set  of  received  pulse  times  is  denoted  by  X  =  {xi,X2,  ■  ■  ■  , xm }  and 
includes  all  detected  and  spurious  pulses.  The  number  of  pulses  in  this  set  is  denoted 
by  |A|  =M. 


2.2  Gaussian  scores 

The  objective  of  this  section  is  to  develop  a  method  for  determining  a  measure  of  consis¬ 
tency  between  a  particular  template  model  (£,  a2)  and  some  observed  pulse  sequence  X. 
Here,  we  begin  by  expressing  the  likelihood  of  making  the  observation  using  the  particular 
Gaussian  model,  in  the  presence  of  the  stationary  channel  noise.  An  exact  expression  for 
such  a  likelihood  function  can  be  expressed  by  enumerating  every  possible  joint  classifica¬ 
tion  of  the  pulses  in  X  in  terms  of  spurious  and  detected  pulses.  In  order  to  facilitate  this,  it 
is  convenient  to  define  the  following  notion: 

An  “association”  (t,x)  G  T  x  X  is  a  2-tuple  that  identifies  a  possible  correspondence  be¬ 
tween  a  template  element  t  G  T  and  an  observed  pulse  x  €  X.  An  “association  relation” 
K  C  T  x  X  is  a  set  of  associations  that  identifies  a  one-to-one  correspondence  between  a 
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subset  of  elements  of  template  T  and  a  subset  of  observed  pulses  in  X.  Each  observed 
pulse  should  be  associated  with  no  more  than  one  template  element,  and  vice  versa.  Some 
template  elements  and/or  pulses  may  not  be  associated  at  all.  The  number  of  associations 
in  the  relation  K  is  denoted  by  |Aj  <min{|X|,|r|}. 


The  likelihood  density  for  observing  X  given  that  “T  starts  at  time  offset  x”  using  the  above 
Gaussian  model  with  stationary  channel  impairments  is  given  by 


p^(X\T,G2,Pd,ps-,x)  =  ]T[  n  -7T^e~i,+:j 

K  \t,x)£K  V-710' 


(2) 


Here,  the  sum  extends  over  every  possible  association  relation  K  between  X  and  T ;  the 
product  kyr_K  extends  over  every  association  pair  within  the  relation  K.  The  likelihood 
function  given  in  (2)  may  be  interpreted  as  an  exact  “score”  with  which  to  measure  the 
consistency  of  the  observation  X  with  the  hypothesis  that  the  template  T  starts  at  time 
offset  x. 


Several  of  the  factors  in  (2)  arc  independent  of  association  relation  K  and  can  be  dismissed 
as  overall  scaling  constants.  After  some  elementary  manipulations,  it  is  found  that  the 
likelihood  density  can  be  decomposed  as  follows: 

p^(X\T,a2,Pd,ps-x)  =  p^(X\T,a2,c0\x)(l  -Pd)lT|piVpsD 


where 

P%{X\T,C2 ,co',l)  =  ^exp(—  c(/,x,t)),  (3) 

K  ( t,x)GK 


and 


c{t,x,x) 


(t  +  x-x)2 

2  a2 


+  CQ  , 


cq  =  In 


Hereafter,  it  is  assumed  that  the  parameters  have  been  chosen  in  such  a  way  so  that  co  <  0; 
this  is  generally  true  unless  the  pulse  environment  is  extremely  dense  or  the  signal  is  highly 
jittered.  The  expression  c(t,x,x)  is  known  below  as  the  “cost”  of  the  association  (t,x) 
for  a  particular-  value  of  the  time  offset  x.  Ah  costs  and  probabilities  appearing  in  (3)  are 
dimensionless. 


It  is  not  practical  to  enumerate  every  association  relation  in  order  to  evaluate  the  above 
likelihood  functions  since  (i)  the  number  of  association  relations  is  exponential  in  both  A 
and  |  T  |  and  (ii)  for  a  given  value  of  the  time  offset  x  most  of  the  association  relations  have 
very  large  values  of  the  total  cost  and  do  not  contribute  significantly  to  (3).  As 

a  result,  a  sub-optimal  strategy  is  needed  to  appropriately  estimate  the  score  as  a  function 
of  x. 
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An  alternative  scoring  function  that  could  be  used  to  approximate  (3)  is  given  by 

q^(X\T,c2,c0;x)  =max  exp(-  ^  c(t,x,x))  .  (4) 

(. t,x)eK 

The  score  given  in  (4)  represents  a  lower  bound  for  (3), 

q^(X\T,o2,c0\x)  <  p^(X\T,a2,c0-,x) . 

The  score  in  (4)  is  known  hereafter  as  the  “q^- score”  and  is  more  often  represented  by  its 
logarithm.  It  may  be  of  interest  to  note  that  the  q^- score  is  restricted  to  the  interval 

0  <  ln^5y;(Z|r,a2,co;x)  <  |T||c0| . 

2.3  Scoring  intervals 

The  q^- score  (4)  has  the  advantage  of  having  a  finite  representation  that  can  be  computed 
without  further  approximation.  To  see  this,  consider  some  particular  scenario  described  by 
(X,T,a2 ,co),  and  define 

K(x)  =  argmin  ^  c(t,x,x ) 

K  ( t,x)€K 

to  be  the  optimal  association  relation  between  X  and  T  at  time  offset  x.  This  relation-valued 
function  is  necessarily  piecewise  constant  in  the  sense  that  there  must  exist  a  partition  of  the 
real  line  M  into  intervals  characterized  by  a  finite  set  of  partition  points  =  { xv  | .v  =  0 ...  5} 
labeled  monotonically 

X0  <Xi  <x2  <  <xs, 

so  that  K(x)  is  unchanged  over  each  interval  of  the  form  x  <G  (x.s_i,x.s).  This  relation  is 
denoted  by  Ks  so  that 

K(x)  =  Ks  for  x  e  (x.s_i ,  x.v) . 

Since  T  and  X  arc  finite  sets,  it  can  be  shown  that  the  partition  can  be  chosen  so  that 
S<2|Z||T|  and 

K{x)  =  0  for  x  <  x0  or  x  >  xs . 


Correspondingly  the  logarithm  of  the  q^- score 


lnq^(X\T,G2,c0-,x) 


-Em6Isc(M4)  if  x  £  (xlS_i,xlS),  s  =  1 . .  .S', 
0  if  x  <  Xo  or  x  >  x5 , 


is  a  quadratic  function  of  x  over  each  interval. 

In  order  to  illustrate  these  concepts  consider  case  in  which  the  observations  and  template 
arc  respectively  given  by 


X  =  {-2,2.1,10,26,31.4,41,69.2,70,73.3}, 

T  =  {0,30,70}.  (5) 


4 


DRDC  Ottawa  TM  2004-220 


The  remaining  parameters  are  taken  to  be  a2  =  25  ps2,  ps  =  0.02/ps  and  PL\  =  0.9,  leading 
to  co  =  —3.581.  The  c/.^'-score  is  illustrated  in  Figure  1(a)  and  the  partition  points  arc 
illustrated  by  vertical  dashed  lines.  The  horizontal  dotted  grid  lines  indicate  multiples  of 
|co|.  (Figures  1(b),  1(c)  and  1(d)  will  be  discussed  in  subsequent  sections.) 

The  segmentation  of  the  real  line  into  intervals  (xv  i .  xiV)  can  be  done  using  a  straightforward 
inductive  algorithm  with  linear  complexity.  The  algorithm  is  discussed  in  in  Annex  A. 

3  Uniform  scoring  profiles 


The  reasoning  described  in  the  preceding  section  may  also  be  applied  to  an  probabilistic 
template  based  on  uniform  distributions.  In  this  case,  we  suppose  the  template  consists 
of  a  pulse  template  T  and  some  time  interval  parameter  w.  As  in  the  previous  section,  T 
is  the  set  of  nominal  times,  but  now  w  is  a  positive  parameter  describing  the  width  of  a 
uniform  distribution  interval.  In  this  case,  instead  of  (1),  the  transmitted  time  variables  are 
distributed  by 

Xj  ~  U(tj  -  w/2  +  x,  tj  +  w/2  +  x) .  (6) 

In  other  words,  xj  has  an  equal  probability  of  taking  on  any  value  that  satisfies  \tj  +  X  —  xj\  <  w/2. 


In  this  case,  we  can  replace  (2)  by  a  likelihood  density  based  on  the  uniform  jitter  model, 
pv(X\T,w,Pd,ps-x)  =  H  w~xl(/t  +  x—x\  <  w/2)]  (7) 

K  (t,x)eK 

x(i-pd)m-WpWpM-iV|>'c 

where l(/\t  +  x  —  x\  <  w/2)  represents  the  indicator  function 

_ / 1  i  f  1  if  If  +  x  —  x\  <  w/2  , 

i(\t+x-x\  <w/2)  ~  1 0  if 


For  simplicity,  it  is  assumed  that  the  elements  of  the  template  T  are  well  separated  compared 
with  the  width  w  so  that,  for  each  value  of  x  —  v  there  is  no  more  than  one  template  element 
t  £  T  that  satisfies  \t  +  x  —  x\  <  w/2. 


Now,  for  x  £  M  and  t  £  T.  let  /,  ( x)  denote  the  number  of  received  pulses  x  €  X  that  satisfy 
\t  +  X  —  x\  <  w/2.  It  is  not  difficult  to  demonstrate  that  (7)  can  be  rewritten  as 

p.u(A|r,w,Fd,ps;x)  =  pu{X\T,w,d0;x){l  -  Pdflp[Xle~p*D 


where 


In  p<u(X\T,w,  d0', x)  =  ^ln 

teT 


PM*)  ) 

(1  -Fd)wps 


(8) 


Here  it  has  been  recognized  that  the  dependence  of  the  right  hand  side  of  the  this  equation 
on  the  parameters  Fd,  ps  and  w  can  be  summarized  by  a  single  composite  parameter,  denoted 


DRDC  Ottawa  TM  2004-220 


5 


X  /  |1S 


(b)  p<ii- score  lnpu(X\T,w,do\x)  using  w  =  5  ,  do  =  —4.511. 


-80  -60  -40  -20  0  20  40  60  80 

x  /  |1S 


(c)  q<u- score  \nqu(X\T.  w,do\%)  using  w  =  5,  do  =  —  4.51 1. 


A(j— 0.5)  /  (is 

(d)  Cross-correlation  histogram  \Rt,xA'j\  usinS  A  =  5 

Figure  1:  Scores  for  the  example  configuration  (5). 
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by  do.  For  later  convenience,  this  parameter  is  taken  to  be 


do  =  —  In 


Pd 

(1  -Pd)wps 


<0. 


Since  lt(x)  is  a  discrete-valued,  and  therefore  piecewise  constant,  function  of  x,  so  must  be 
In pu.  That  is,  it  must  be  possible  to  partition  the  real  line  at  a  finite  number  of  points  so 
that  lnpu(X\T,w, do'.x)  is  constant  over  each  interval.  In  fact,  the  set  of  discontinuities  in 
In  pu  is  almost  surely  equal  to  the  finite  set 


‘Pu  =  {x  £  K|  3x  £  X,  t  £  T,  x  =  v  —  t  ±  w}  . 


Consequently,  a  brute-force  approach  to  computing  the  complete  partition  is  to  first  calcu¬ 
late  each  of  the  2\X\ |T|  differences  x  —  t±  vv,  and  then  sort  the  resulting  set. 

Figure  1(b)  illustrates  the  the  pu- score  for  the  example  in  (5)  using  parameter  values  vv  =  5 , 
ps  =  0.02,  F\i  =  0.9.  Note  that  the  dotted  horizontal  grid  lines  represent  multiples  of  | z/q |  = 
4.511. 


If,  instead  of  the  summation  YjK  'n  (7),  a  max^  operation  is  used,  then  the  result  would 
become 


qu(X\T,w,Pd,ps',i)  =  max  TT  [w  ll(\t  +  x -x\  <  w/2)\ 

{tjc)eK 

xll-P/H^pF'-'V0 

=  qv(X\T,w,do;z)0 -Pa)|I'lpiX|c_l’,D. 


where 

\nqv(X\T,w,do\x)=N(x)\do\.  (9) 

Here  N(x)  indicates  the  number  of  template  elements  t  G  T  for  which  there  is  at  least  one 
pulse  x  £  X  satisfying  \t  +  x  —  x\  <w  /  2.  Like  In  pv,  this  function  is  piecewise  constant;  its 
discontinuities  form  a  subset  of  Tw. 

The  score  in  (9)  is  a  lower  bound  for  (8) 

qu(X\T,w,d0-,x)  <  pu(X\T,w,d0-,x) , 

and  In q ‘u  is  restricted  to  the  interval 

0  <  lnqu(X\T,a2,c0\x)  <  \T\\do\  ■ 


The  q.y-scorc  is  illustrated  in  Figure  1(c)  for  the  example  (5)  using  parameter  values  vv  =  5 , 
ps  =  0.02,  Po  =  0.9  and  do  =  —4.511. 
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4  Histogram  scoring  profiles 


In  Section  5,  the  q.^-  and  py-scoring  techniques  arc  contrasted  against  the  cross-correlation 
histogram  presented  in  [5].  The  latter  technique  is  briefly  paraphrased  here.  In  that  ap¬ 
proach,  the  received  pulse  time  set  X  is  compared  to  the  template  set  T  using  the  cross¬ 
correlation  function 

x€X  t<=XT 

where  8(-)  refers  to  the  delta  function  of  Dirac.  This  function  is  then  smoothed  in  some 
fashion.  In  [5],  the  smoothing  is  accomplished  by  integrating  over  sequential  intervals  of 
the  form  ((j  —  1)A,  jA]  where  A  is  some  specified  bin  width  and  j  is  a  bin  counter.  The 
result  is 

/\/'A 

/  Srjci't)  dx  =  \Rt,x,aj\  > 

J{J-  1)A 

where  \Rt.x,A;j\  counts  the  number  of  elements  in  the  set 

Rt,x,aj  =  {{t,x)\t  G  T,xeX,(j-  1)A  <x-t  <  jA}  . 

Like  the  association  relations  introduced  in  Subsection  2.2,  each  set  Rt.x,aj  C  T  x  X  is  a 
collection  of  associations  between  elements  of  T  and  X.  However,  unlike  the  previous  rela¬ 
tions,  there  is  no  one-to-one  requirement  on  Rt,x,aj  so  that  two  or  more  different  pulses  in 
X  may  be  associated  with  the  same  template  element  in  T.  Although  the  perfect  association 
would  have  exactly  |T|  elements,  there  is  no  fundamental  upper  bound  on  \Rt.x.a-.j\- 

By  plotting  |/?r,x,A;./i  against  the  bin  centers  (j  —  0.5 )A,  for  sequential  values  of  j,  one 
obtains  the  cross-correlation  histogram  of  described  in  [5]. 

For  illustration,  the  cross-correlation  histogram  for  the  configuration  in  (5)  is  shown  in 
Figure  1(d)  using  a  bin- width  of  A  =  5.  The  optimal  number  of  associations  is  |T|  =  3  and 
is  indicated  by  a  solid  horizontal  line.  Note  that  the  bin  corresponding  to  x  G  (0,5]  contains 
|^r,x,A;t|  =  4  associations,  despite  T  containing  only  three  elements.  This  function  can 
be  compared  with  the  other  scoring  functions  shown  in  Figure  1  for  the  same  example 
configuration. 


5  Simulation 


This  section  demonstrates  the  effectiveness  of  the  scoring  approaches  introduced  here  and 
compares  them  with  the  cross-correlation  histogram  of  [5]  by  using  a  synthetic  data  set. 

5.1  Configuration 

Following  the  example  of  [5],  a  synthetic  experiment  is  conducted  in  which  three  signals 
are  interleaved  and  observed  through  a  noisy  channel. 

The  characteristics  of  these  signals  are  as  follows: 
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1 .  The  first  signal  is  a  set  of  transmitted  pulses  X\  described  by  a  constant  pulse-to-pulse 
interval  of  875ps,  without  jitter. 

2.  The  second  signal  X2  is  a  5  element/5  position  staggered  sequence  whose  pulse-to- 
pulse  intervals  arc  { 620/,/s.  325/,/s.  220/rs.  475/./S.  490/,/s}.  These  intervals  arc  repeated  in 
cyclic  order,  without  jitter 

3.  The  third  signal  X2  exhibits  non-cumulative  jitter  with  a  nominal  pulse-to-pulse  interval 
of  875/,/s.  The  actual  transmission  time  of  each  pulse  is  a  random  variable  normally 
distributed  about  a  nominal  time  with  a  standard  deviation  of  1 2ps. 

In  the  simulation,  each  of  these  signals  is  persistent  over  a  long  period  of  observation  and 
each  is  characterized  by  a  starting  phase  that  has  been  selected  from  a  uniform  distribution. 

The  complete  set  of  transmitted  pulses  is  the  union  X  =  X1UX2UX3.  This  sequence  is  then 
mapped  into  a  received  pulse  sequence  X  1— >  X  by  a  stationary  noise  model  in  which  (i) 
each  pulse  in  X  has  a  probability  of  0.9  of  being  detected  and  appealing  in  X,  and  (ii)  X 
includes  additional  spurious  pulses  that  arc  distributed  according  to  a  Poisson  process  over 
the  observation  width  with  average  density  of  1  pulse  every  1 5 Ops.  The  received  set  X  is 
the  union  of  all  detected  and  spurious  pulses. 

5.2  Results 

Several  scores  for  this  scenario  arc  computed  for  values  of  the  offset  in  the  interval  0  <  x  < 
6000,i/s.  These  arc  illustrated  in  Figures  2-7.  In  each  figure,  the  three  panels  show  q^  ~,  pu~ 
and  histogram-scores  with  the  horizontal  dotted  grid  lines  indicating  multiples  of  |co|,  |do| 
or  1  respectively.  The  i/?/-scorcs  have  been  omitted  since  they  arc  almost  indistinguishable 
from  the  //^-scores.  Where  relevant,  the  scoring  parameters  ps  =  0.02/is  and  P&  =  0.9 
have  been  used  to  compute  c 0  and  do. 

Figures  2  and  3  show  the  results  of  comparing  X  with  a  pulse  template  containing  eleven 
equally  spaced  elements, 

7j  =  {0.0,875.0,1750.0,2625.0,3500.0,4375.0, 
5250.0,6125.0,7000.0,7875.0,8750.0}, 

corresponding  to  the  expected  behaviour  of  signal  X\ .  In  Figure  2,  the  q^-  and  pv- scores 
arc  calculated,  respectively,  using  a  =  3,i/s  and  w  =  3,i/s.  The  histogram  is  computed  using 
A  =  3ps.  Figure  3  illustrates  the  q^-,  pw  and  histogram-scores  using,  respectively,  a  = 
1 2,t/s,  w  =  1 2,t/s  and  A  =  1 2,i/s.  It  should  be  noted  that,  since  there  is  no  jitter  on  the  signal 
X\ ,  the  scoring  approaches  with  the  smaller  values  of  a,  w  or  A  are  expected  to  be  more 
suitable  than  those  with  the  larger  values. 

Figures  4  and  5  show  similar  results  using  the  template 

T2  =  {0.0,620.0,945.0,1165.0,1640.0,2130.0, 

2750.0, 3075.0, 3295.0, 3770.0, 4260.0} , 
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x/ps 

(b)  p^-score  with  w  =  3ps  ,  do  =  —5.017. 


A(j— 0.5)  /  |1S 

(c)  Histogram  with  A  =  3. 

Figure  2:  Cross-correlation  scores  forX  computed  with  template  T\  and  narrow  tolerance. 
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T/|1S 


(b)  p<u- score  with  w  =  12/rs  ,  do  =  —3.651. 


A(j— 0.5)  /  (is 


(c)  Histogram  with  A  =  12. 

Figure  3:  Cross-correlation  scores  forX  computed  with  template  T\  and  wide  tolerance. 
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corresponding  to  the  expected  behaviour  of  signal  X^.  Like  the  previous  signal,  X3  has  no 
jitter  on  the  arrival  times,  and  the  smaller  values  of  a,  w  or  A  shown  in  Figure  4  should  be 
considered  to  more  suitable  than  the  larger  values  shown  in  Figure  5. 

Finally,  Figures  6  and  7  show  the  scores  when  the  template 

T3  =  {0.0,600.0,1200.0,1800.0,2400.0,3000.0, 

3600.0, 4200.0, 4800.0, 5400.0, 6000.0}  , 

is  applied.  This  template  represents  the  expected  behaviour  of  X3.  Since  the  signal  is 
associated  with  a  non-cumulative  la  jitter  of  1 2/,/s.  the  smaller  values  a  =  3//s,  w  =  3/,/s  or 
A  =  3/,/s  shown  in  Figure  6  are  not  expected  to  perform  as  well  as  the  larger  values  a  =  1 2/,/s, 
w  =  1 2p s  or  A  =  1 2p s  shown  in  Figure  7. 

Some  observations  about  these  results  are  in  order.  Both  X\  and  X3  represent  idealized 
signals  that  contain  no  jitter,  or  time-of-arrival  uncertainty.  Although  this  situation  is  some¬ 
what  idealized,  it  illustrates  the  limiting  behaviour  of  the  various  scoring  techniques.  When 
the  templates  7j  and  73  are  applied  to  the  interleaved  signal,  we  see  that,  in  all  cases,  it  is 
preferable  to  use  smaller  values  of  the  parameters  a,  w  and  A.  Moreover,  the  histogram 
performs  as  well  or  better  than  the  q-  and  qu  scores. 

The  signal  X3  is  more  realistic.  It  is  fairly  clear  that  the  larger  parameter  shown  in  Figure  7 
result  in  a  cleaner  signal  detection  than  do  the  smaller  values  shown  in  Figures  6.  In  fact, 
the  A  =  3qs  histogram  results  in  the  most  ambiguous  results.  By  contrast,  the  cleanest  and 
most  consistent  results  are  generated  by  the  q^- score  with  a  =  1 2/,/s.  possibly  owing  to  the 
fact  that  this  score  is  a  most  closely  matches  the  signal  generation  model. 


6  Conclusions 


This  paper  has  described  candidate  techniques  that  may  be  applied  in  Electronic  Support 
for  identifying  fixed  patterns  of  pulses  in  an  observed  sequence  specified  using  an  a  pri¬ 
ori  template.  The  techniques  described  apply  to  temporal  pulse  sequences  in  which  the 
pulse  descriptor  words  consist  of  pulse  arrival  times.  However,  it  is  not  difficult  to  modify 
these  techniques  to  include  more  general  pulse  parameters,  such  as  carrier  frequency,  and 
pulse  modulation  parameters  whose  measured  values  may  follow  some  known  probability 
distribution. 

The  first  method,  known  herein  as  the  q^- score,  attempts  to  formulate  a  lower  bound  for 
the  likelihood  of  observing  a  measured  sequence  given  some  a  priori  model.  The  com¬ 
plete  model  consists  of  a  characteristic  template  of  relative  pulse  arrival  times,  a  Gaussian 
non-cumulative  jitter  variance  and  two  stationary  parameters  describing  the  spurious  pulse 
density  and  pulse  detection  probability.  An  algorithm  is  presented  whereby  an  exact  repre¬ 
sentation  for  the  r/.y-scorc  can  be  evaluated  with  linear  time  complexity. 

The  second  method,  known  as  /;?/-score,  is  similar  to  the  r/.y-score,  but  attempts  to  compute 
the  likelihood  of  observing  a  sequence  using  a  non-cumulative  jitter  model  that  follows  a 
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T  /  (J.S 


A(j— 0.5)  /  |IS 

(c)  Histogram  with  A  =  3. 

Figure  4:  Cross-correlation  scores  forX  computed  with  template  T2  and  narrow  tolerance. 
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"t/gS 


(b)  p<u- score  with  w  =  12/rs  ,  do  =  —3.651. 


A(j-0.5)  /  (is 

(c)  Histogram  with  A  =  12. 

Figure  5:  Cross-correlation  scores  forX  computed  with  template  T2  and  wide  tolerance. 
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(b)  pzi- score  with  w  =  3/rs  ,  do  =  —5.017. 


2.5 


A(j— 0.5)  /  |IS 

(c)  Histogram  with  A  =  3. 


Figure  6:  Cross-correlation  scores  forX  computed  with  template  r3  and  narrow  tolerance. 


DRDC  Ottawa  TM  2004-220 


15 


(b)  p<u- score  with  w  =  12/rs  ,  do  =  —3.651. 


6 


A(j— 0.5)  /  |ts 

(c)  Histogram  with  A  =  12. 


Figure  7:  Cross-correlation  scores  forX  computed  with  template  r3  and  wide  tolerance. 
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uniform  distribution.  A  lower  bound  for  the  pu- score,  known  as  the  r/.y- score,  is  also 
presented  and  is  found  to  be  nearly  identical  to  the  pu- score  in  practice. 

Each  of  the  q^-  and  py-scores  depend  on  two  parameters.  One  parameter  (either  a  or  w 
respectively)  measures  the  tolerance  of  the  arrival  time  around  some  expected  value.  The 
second  parameter  (either  co  or  do)  measures  the  cost  of  associating  elements  of  the  template 
with  pulses  in  the  observation  and  depends  on  the  estimated  pulse  density  and  detection 
probability. 

The  two  methods  arc  compared  against  a  cross-correlation  histogram  introduced  in  [5].  It 
is  found,  using  a  synthetic  experiment,  that  for  signals  without  jitter  or  time-of-arrival  un¬ 
certainties,  the  ry^'-and  py-scorcs  have  similar  or  slightly  worse  discriminating  power  than 
the  cross-correlation  histogram.  However,  for  a  signal  with  a  significant  degree  of  jitter, 
the  <y jy  - score  whose  parameters  are  chosen  to  match  the  expected  jitter  variance  provides 
a  cleaner  template  detection  measure  than  either  the  py-score  or  the  cross-correlation  his¬ 
togram. 
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Annex  A 

Segmentation  Algorithm 


In  Subsection  2.3,  it  was  claimed  that  the  g^-score  defined  in  (4)  is  piecewise  smooth  and 
quadratic  over  intervals  of  the  form  (x.s_i  ,x.v)  where  ‘L^r  =  {xs|.v  =  0. .  .5}  represents  some 
collection  of  partition  points.  This  appendix  describes  an  efficient  recursive  algorithm  with 
which  one  can  determine  the  elements  of  T^,  in  increasing  order. 

Suppose  that  some  partition  point  xs_i  is  known  and  it  is  required  to  find  the  next  subse¬ 
quent  partition  point  ts.  Within  the  interval  x  e  (kv-i,Tv),  the  optimal  association  relation 
Ks  is  assumed  to  have  been  determined.  For  each  t  £  T,  there  arc  two  possibilities: 

1.  It  may  be  that  for  all  x  £  X,  the  cost  of  associating  t  with  x  is  non-negative  c(t,x,x)  >  0 
for  Vx  £  (xv_  i  ,xv).  If  this  is  the  case,  then  the  optimal  association  relation  Ks  does  not 
include  any  association  for  t  and  we  may  say  that  t  is  “unassociated”  on  this  interval. 

2.  Alternatively,  there  may  exist  x  £  X  such  that  c(t,x,x)  <0forx£  (xv  |.xv).  In  this  case, 
we  can  define  x  =  argminrGXc(t,x,x).  Then,  (t,x)  £  Ks  and  t  is  said  to  be  “associated 
withx”  on  (x^-ijXs). 

Now,  for  a  particular  value  of  t  £  T,  it  is  of  interest  to  find  the  smallest  value  of  x  exceeding 
x.s_i  at  which  the  association  for  t  within  K(x)  exhibits  a  change.  This  is  denoted  by  xs_t 
and  is  guaranteed  to  satisfy  xs_i  <  xv  <  XSJ.  The  criterion  for  determining  xA.,  depends  on 
whether  or  not  t  is  associated  in  Ks. 


1.  If  t  is  unassociated  on  (xv  i.xv),  then  xS)f  is  the  first  value  beyond  xs_ i  at  which  it 
becomes  beneficial  to  add  an  association  for  t.  More  precisely,  we  have 


X, ,  =  inf< 


|x  >  Xi_i  3x,c(t,x,x)  <  o| . 


2.  However,  if  t  is  associated  with  some  x  on  (x.v .  i  .xv),  then  xsA  is  the  first  value  after 
xv  i  at  which  it  becomes  beneficial  either  to  remove  the  association  (t,x)  or  to  replace 
it  with  a  different  association  for  t, 

zsA  =  infjx  >  x.v_i  | Vx,c(t,x,x)  >  o|  U  |x  >  Xj_i \3x,c(t,x,x)  <  c(t,x,x)| . 

For  sufficiently  large  values  of  XiS_  i ,  some  values  of  zsj  may  not  exist  since  the  relevant 
infimum  may  extend  over  an  empty  set.  If  none  of  these  values  exists,  then  we  may  conclude 
that  there  are  no  partition  points  exceeding  xiV-i  so  that  5  =  5—1.  However,  if  some  of  the 
infima  exist  then 


Is  =  rnin{xS),} 

t  ’ 

ts  =  argmin  {xu  } 

t 
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It  is  then  straightforward  to  determine  from  Ks,  depending  on  which  of  the  above  cases 
gave  rise  to  the  relevant  infimum,  by  either  adding,  removing  or  replacing  the  association 
involving  the  ts. 

For  the  implementation  of  the  algorithm  discussed  below,  it  is  required  to  assume  that  the 
elements  of  T  =  {tn\n  =  1 . . . N}  and  X  =  {xm\ m  =  1 . . . M )  arc  indexed  in  increasing  order 
so  that  order  so  that  t\  <  t2  <  ■  ■  ■  <  W  and  x\  <  X2  <  ■  ■  ■  <  xm-  As  discussed  above,  the 
algorithm  operates  on  the  premise  that  a  particular  tn€T  need  be  associated  only  if  there 
exists  some  r£lso  that  c(t.x.z)  <  0.  Referring  to  (3),  it  is  found  that  this  condition  is 
equivalent  to 

t~  <x-x  <f+ 

where  t+  =  tn  ±  \/ —2a2 c\).  It  is  furthermore  assumed  that  these  thresholds  do  not  overlap, 
ensuring  that 

?r  <t\  <f+  <  q  <t2<  t+  <  ?3  <  . . .  <  t+_j  <  t~  <  tn  <  t+ . 

When  applied  methodically,  these  variables  can  be  used  to  evaluate  the  infima  described 
above. 

For  convenience,  the  algorithm  has  been  divided  into  several  procedures  shown  in  Fig¬ 
ures  A.  1,  A.2,  and  A.3.  The  entry  point  is  the  method  “IteratePartition”  that  takes,  as 
parameters,  all  values  of  tn,  and  xm.  Although  the  only  explicit  output  of  this  algorithm 
is  the  set  of  partition  points  '2^,  additional  processing  can  be  applied  on  an  interval-by- 
interval  basis  by  inserting  an  optional  method  called  “ProcessAssociationRelation”. 


DRDC  Ottawa  TM  2004-220 


19 


procedure  =  IteratcPartition( { (t~  . tn , t^)\n  =  1 . .  .N},{xm\m  =  1 . . . M }) 

for  Vn  =  1 ...  A  do 

4—  1 

qn  4—  Unassociated 

4 

end 

N^N 
n^N 
<1  -  {x,t} 

repeat 

if  q„  =  Associatedlmprove 
<—  m„  +  1 
Make  Association 
elseif  qa  =  Unassociated 
Make  Association 
elseif  q,',  4—  AssociatedLose 
RemoveAssociation 
if  A  =  0  return  end 
end 

h  4—  argminj;{x„|n  =  1 . .  .A} 

^A  U  {X/;} 

ProcessAssociationRelation 

end 

end 

Figure  A.1:  “IteratePartition”  iterates  through  the  interval  partitions  T^-  =  { xA  |  >y  =  0. .  ..S’}, 
procedure  MakeAssociation 

if  mu  <  M  do 

t  4  (Xthjj  hr  Xihf,  i  I  )  /  2  tft 

if  x  <  x  do 

qn  =  Associatedlmprove 
t 

return 

end 

end 

qn  =  AssociatedLose 
X/5  —  X_ 

return 

end 

Figure  A.2:  “MakeAssociation”  determines  the  state  parameters  for  the  current  interval 
after  an  association  has  been  made. 
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procedure  RemoveAssociation 
q„  =  Unassociated 
inn  <—  mh  +  1 
if  rhfi  <  M  do 

Xa  <—  x,n-  — 

else 

end 

return 

end 

Figure  A.3:  “ RemoveAssociation ”  determines  the  state  parameters  for  the  current  interval 
after  an  association  has  been  eliminated. 
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